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ABSTRACT 

The kinetics of monoatomic steps in diffusion-controlled crystal growth and 
evaporation processes are investigated analytically using a Green's function ap- 
proach. Integro-differential equations of motion for the steps are derived; and a 
systematic linear stability analysis is carried out treating simultaneously pertur- 
bations both along and perpendicular to the steps. Morphological fluctuations 
of steadily moving steps in response to ambient thermodynamic noises are also 
studied within a general Langevin formalism. Finally, a phase field model is 
developed to investigate the time-dependent, collective motion of steps. An ap- 
plication of the model to a finite step train recovers a variety of kinetic behaviors 
such as the bunching and spreading of steps. 



PACS numbers: 61.50.Cj, 05.40.+j, 68.35.Ja 



1. INTRODUCTION 



The kinetics of atomic steps play a central role in the growth of singular 
crystal surfaces from vapor, liquid, and solid phases, in a wide variety of ma- 
terials. Step motion is also a predominant elementary process governing such 
macroscopic changes of crystal surface topography as dissolution, faceting, etch- 
ing and coarsening. Numerous investigations, both theoretical and experimental, 
have been devoted to understand the morphology and dynamics of atomic steps. 
Recent resurgence of interests in step dynamics is largely stimulated by the de- 
velopment of crystal growth methods with atomic-precision, and by the advances 
in high- resolution microscopies. 

This paper addresses several theoretical aspects of the kinetics of diffusively 
coupled steps. Our investigation is based upon the classical model of step ki- 
netics developed by Burton, Cabrera, and Frank (BCF)!^' The ideal BCF model 
visualizes the singular crystal surface as consisting of closely-packed terraces sep- 
arated by elementary steps. Growth of crystal then proceeds through the lateral 
progression of these steps by incorporating diffusing adatoms from the terraces. 
The BCF model, supplemented with later modifications and extensions, has been 
extensively investigated and remains the most mathematically well-posed model 
of crystal growth by step flow dynamics. Although we restrict ourselves to the 
behavior of monoatomic steps governed by surface diffusion, and to situations 
where other assumptions of the BCF model are valid, we hope that this study 
can stimulate progress in elucidating the role of steps (not necessarily elementary) 
under other growth conditions, e.g., the growth of crystal by molecular epitaxy, 
or the ledge growth by volume diffusion in solid/solid phase transformations. 

One central question in step kinetics concerns the morphological stability of 
step patterns. Steps are seldom seen to be equally spaced. It is also observed 
that fluctuations in the distance between steps can often lead to the bunching 
of them into multi-step bands. Furthermore, meandering of the steps along the 
step direction on scales much greater than molecular lengthscale is common- 
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place. Bales and Zangwill recently pointed out that, in addition to the effect 
of equilibrium fluctuations, a diffusive instability can contribute to the mean- 
dering of steps under growth situations. To address the stability issue, linear 
stability analyses, treating separately fluctuations in the spacing between adja- 
cent steps and that along the step direction, have been carried out. Nevertheless, 
two popular approximations are often adopted. First, when considering pertur- 
bations in the step spacings, it is commonly assumed that the velocity of a step 
only depends upon the width of its two adjacent terraces. Second, in treating 
fluctuations along the steps, and in studying the motion of steps in general, a 
quasi-static approximation is used by assuming sufficient slow motion of steps. 
While these approximations are believed to be good for widely separated steps 
and under small driving forces, a fully time-dependent analysis is warranted in 
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the opposite situations. As emphasized by Ghez and Iyer, such situations may 
arise particularly in certain epitaxy experiments where the so-called fast steps are 
encountered and where nearby steps are strongly coupled by their overlapping 
diffusion fields. Linear stability analysis which takes account of the overlap of dif- 
fusion fields has recently been carried out numerically!*"'' In this paper, the step 
fiow problem is formulated using a Green's function approach. Time-dependent 
integro-differential equations of motion for the step positions are derived, from 
which a systematic linear stability analysis is carried out. This formulation al- 
lows a simultaneous treatment of generic perturbations both along step extension 
and in the step spacings. The validity of the quasi-static approximation and the 
effect of the overlap of terrace diffusion fields are also examined. These works are 
discussed in section 2 through section 4. Section 5 analyzes the effect of asym- 
metric attachment kinetics on the linear stability of infinite step trains, under 
both growth and evaporation situations. 

In section 6, we discuss the morphological fiuctuations of steps during non- 
equilibrium growth using a Langevin formalism. In this formalism, local ther- 
modynamical noises on the terraces and on the steps are assumed to be the sole 
sources of fluctuations which get amplifled dynamically during growth and lead 
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to the roughness of steps. We show that this method consistently recovers the 
equihbrium fluctuations. Whilst step roughness in equilibrium does not depend 
on step spacing, we found that the diffusive coupling of steps in motion during 
growth gives rise to a correlation between fluctuations on different steps. 

While stability analysis gives important information with regard to the role 
of step kinetics in crystal growth, more insights can only be provided by detailed 
investigations of the temporal development of step configurations. It is crucially 
required, for example, to interpret in-situ experimental observations. In section 
7, we develop a phase field approach to study the time-dependent motion of finite 
and generically non-equidistant step trains. This method is numerically efficient, 
easily adaptable to two-dimensions, compatible with shape fluctuations along the 
step extension, and able to simulate complicated time-dependent morphology of 
step patterns. 



2. THE BCF MODEL 



Consider the mathematical formulation of the original BCF model, for a 
close-packed crystal surface below its roughening transition temperature, grow- 
ing from a supersaturated mother medium. We concentrate on the layer-by- 
layer growth mechanism characterized by the diffusion-controlled propagation of 
monoatomic steps, and neglect entirely the effect of two-dimensional island nu- 
cleation, under the assumption of small supersaturations. For concreteness, it 
is further assumed that the growth units are incorporated to the steps through 
diffusion along the terraces only. Although this mode of building growth units 
is predominant in vapor and solution growth, the alternative mode of depositing 
growth units directly from the bulk mother medium is important in certain solu- 
tion growth environments and particularly in solid-solid phase transformations. 
The BCF formulation is equally applicable to crystal evaporation and dissolution 
processes. 
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In the BCF model, transport of adatoms on the terraces proceeds through 
surface diffusion according to 




c — 



(2.1) 



where c(r, t) is the adatom density, is the mean hfetime of adatom on the ter- 
races before evaporating into the vapor, c^o is the adatom density far away from 
the steps, and D is a diffusion coefficient assumed isotropic over all terraces. In 
the present macroscopic description, steps are regarded as mathematically sharp 
dividing lines between terraces. A crucial assumption, inherent in most contin- 
uum treatments, demands that the monoatomic steps are molecularly rough with 
numerous kinks excitations so as to behave as ideal line sinks for growth units. 
The condition of molecular roughness of steps is also a prerequisite for us to safely 
neglect lateral diffusion of adatoms along the steps, and to assume isotropic ad- 
vancement of steps. A lower kink density hence smoother steps often implies 
stronger effect of attachement kinetics and more pronounced anisotropy in step 
motion and morphology. These and other assumptions of the BCF approach 
must be carefully examined when applying the original model to experiments, 
particularly to such nonconventional situations as molecular beam epitaxy'"^' at 
low temperatures. 

Relevant microscopic details of the atomic processes near the steps are re- 
flected in the boundary conditions on the sharp steps. Firstly, conservation of 
materials requires that the normal velocity of a step satisfies 



where Cp is the change of atomic density in the area swept by the moving step, 
n is the unit normal to the step, and (Vc)± are gradients computed on the two 
sides of a step, respectively. Secondly, another boundary condition relates the 
density of adatoms at the step to its local equilibrium value, taking into account 



CrVn = D [(Vc)+ - (Vc)_] • n 



(2.2) 
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the curvature correction via the Gibbs-Thompson relation: 



Cstep = Co I 1 + 



7O 




(2.3) 



where cq is the adatom density for a straight step in equiUbrium, 7 is the isotropic 
hne tension of the step, Q is the atomic area of the sohd, and k is the step 
curvature. The density of adatoms far away from the steps is maintained above 
the equihbrium density at Cqo > cq through supersaturation, providing the driving 
force for step propagation. We have neglected the effect of crystalline anisotropy 
and also temporarily ignored the elevation of adatom density at the steps due 
to attachment kinetics. Equations (2.1) through (2.3) complete the description 
of the symmetric BCF model for crystal growth by the step flow mechanism, 
where atomic processes on both sides of a step contribute cquivalently to step 
progression. The corresponding asymmetric version of the model, accounting for 
nonequivalent attachment kinetics as first discussed by Schwoebel and Shipseyj''' 
will be considered in section 5. 

Multiple lengthscales are involved in the current problem, whose relative 
ratio characterize different growth regimes and determine the validity of vari- 
ous approximations. With the assumption of mathematically sharp steps, the 
continuum model is understood to apply on length scales much longer than the 
microscopic capillary length do (to be defined below). Among the relevant macro- 
scopic lengthscales, the mean adatom diffusion length on the terraces is given by 
Xg = \/Drsi ciiid the typical step spacing is denoted by A. Another quantity 
I = 2D/v with dimensionality of length can be conveniently introduced as an 
alternative measure of the speed of step motion. Defining a quantity 



and measuring time in units of Ts, length in units of a^g, and velocity in units of 



u{y,t) 



c{r,t) 



(2.4) 
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D/tsi the equations of motion can be written in the dimensionless form 




(2.5) 



[(Vu)^ - (V^i)-] • n 



(2.6) 



Ustep = A - d^K 



(2.7) 



where A = (cq© — CQ)/cr is the dimensionless supersaturation and cLq = 
^ilco/{kBTcrXs) is the dimensionless capillary length. The problem defined by 
(2.5) through (2.7) resembles closely the problem of solidification, in spite of the 
extra homogeneous term —u on the r.h.s. of (2.5). 



The moving boundary problem (2.5) to (2.7) belong to a general category of 
problems known as the Stefan problem, and is strongly non-linear in character. 
One must solve a diffusion-like equation subject to boundary conditions which 
themselves depend upon the solution of the same diffusion equation. Complete 
analytic treatment of the entire problem in 2-space and 1-time dimensions is 
difficult. Even the linear stability analysis has only been carried out for syn- 
chronously moving equidistant steps, under quasi-static approximations. Studies 
of time-dependent step motion, particularly for non-equidistant steps, are even 
more limited. Fortunately, progress can be made analytically by reformulating 
the full 2-1-1 dimensional problem in terms of a closed set of integro-differential 
equations for the positions of the steps. This boundary integral method, using 
Green's functions, was first put forward by Langer and Turski in their study of 
directional solidification fronts. A straightforward application of their procedure 
to our present problem yields the dynamical evolution equations for the moving 
steps, starting from which many detailed analyses can be made. 



3. EQUATION OF MOTION FOR THE STEPS 
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Consider an infinite train of non-overlapping steps on the crystal surface 
shown schematically in Fig.l. The time-dependent position of the nth step is 
parameterized by z = ('^{x,t). In a coordinate frame moving with velocity V in 
the positive 2;-direction, the free boundary problem takes the form: 

V' + V^^-l-^y{x,z,t) = (3.1) 

(v + dx^- [{Vu)l - {S/u)!] ■ do-" (3.2) 

u{x,C,t) = A-doK{C} (3.3) 

where the superscripts label the steps and the vector step length-element dcr 
is defined according to Fig.l. To proceed, we define a Green's function in the 
moving frame by 

with the short notation p = {x,z,t), pi — {xi,zi,ti). The causality condition 
that G{p\pi) = for t < ti is also imposed. For an arbitrary space-time point p, 
we have from (3.1) and (3.4), 



u{p) ^ J dtij dai ■ [u{pi)ViG{p\pi) - G{p\pi)Viu{pi)] 

- f (3.5) 
+ V j dti j dxidzi-^[u{pi)G{p\pi)] 



-oo Ai 



where the Green's theorem has been applied to the region Ai containing the 
point p but excluding all the steps bounded by the boundary surface ^i. The 
volume contribution in (3.5) vanishes identically because the integration over dzi 
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gives cancelling terms from two sides of each step. This cancellation is specific to 
the symmetric model, where the continuity of adatom density across each step 
is maintained. Further application of this continuity property to the first part of 
the surface integral in (3.5) also gives a null result. Hence (3.5) reduces to 

+00 t „ 

<P)= Yl G{p\pi)[{Viu{pi))--{Viu{pi))+]-dcTi+ (3.6) 

which, upon using the boundary condition (3.2) and letting p approach the nth 
step, leads to the equation of motion for the nth step: 

+00 * 

A-do^{C}= Yl I '^^^ I dxiG{x,C.t\^iXTM){v +'^Y (3.7) 



-00 —00 



Here the curvature is by definition 



2 An 



1 + 



dC 

dx 



-3/2 



(3.8) 



This set of closed, nonlinear, integro-differential equations describe the time- 
dependent evolution of an infinite train of steps mutually coupled through the 
overlapping adatom diffusion fields. With little modification, the boundary in- 
tegral approach can also be used to investigate pattern selection problems for 
spiral steps. Finally, we point out that, in the present consideration, no energetic 
interaction between steps, such as that due to elasticity, is yet incorporated. 
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4. LINEAR STABILITY ANALYSIS 



The boundary integral formulation makes possible certain analytical treat- 
ments otherwise infeasible using the original full diffusion model. For example, 
this approach is particularly suitable for determining various steady state solu- 
tions and for analyzing their stabilities. Although equation (3.7) supports in 
general steady state solutions of both straight and curved steps, we restrict our 
consideration in this paper to straight steps only. 

This section presents a linear stability analysis of uniformly moving step 
trains. The advantages of using (3.7) for linear stability analysis, rather than 
using the full diffusion equation, are two-fold. First the method is systematic and 
allows us to treat simultaneously perturbations both along and perpendicular to 
the steps. Second, the calculation can be carried out without having to use the 
quasi-static approximation, which becomes unreliable for both fast-moving and 
densely-packed step trains. The situation of fast step trains is likely to occur in 
MBE growth where the dimensionless supersaturation is typically large. Some of 
our results reported here were previously known but were obtained using different 
methods under quasi-static approximations. 

First we concentrate on steady state solutions for straight steps. Consider 
an infinite train of equidistant steps moving uniformly with velocity V in the 
positive 2;-direction. In the co-moving frame of the steps, letting ('^{x,t) = nX 
with A being the step spacing, and using the real-space representation of the 
Green's function 



{x-xi)^ + {z-zi+V{t-ti)y 
4(t-ti) 



(4.1) 



we obtain from (3.7), 



+00 +^ 

A = V Y. 



m=—oo 



+00 +00 

J dr J (ia;iG(0,0,T|a;i,mA,0), (4.2) 

-oo 
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which leads to the transcendental equation 



2A 



a-\- — a.- 



where we have defined for convenience 



, Ct-I-A , OL-\ 

coth — h coth — — 



(4.3) 



a±^yi + i±l. (4.4) 

Equation (4.3) can be solved numerically for the velocity of steps V = V^(A, A), 
which is monotonically increasing in both A and A. 

Next let us analyze the morphological stability of this system of equidistant 
steps. Consider small-amplitude perturbations to the steady state step profiles 
in the general form 

C{x,t) = C{.^.t)-n\ (4.5) 

for all n. Equations governing the time evolution of perturbations are derived 
by expanding (3.7) in powers of the small quantity CJ^ . We shall use the Green's 
function in a Fourier representation 

^^'^'^ i {2ny 2M{k,u) ^ ^ 

where 



M{k,uu) = ^k^ + 1 + ^ + 100 (4.7) 

with a positive real part. The integrand in the Green's function can be readily 
expanded to first order in perturbation to read 

where the zeroth order term is given by 

= g-(n-m)A//-|n-m|AM (49) 
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and the first order term takes the form 



-{C-(i)/l-\C-Ci\M forn = m 

- {Ml + l)e-(^+l/0(n~m)A^^n _ ^^my^ forn > m 



(4.10) 



Now these expressions can then be inserted into the r.h.s. of (3.7) which has the 
expansion 



/ 



dtidxi ,^ ^r, y. 

m 



{2ny 



2M 



dti 



(4.11) 



It is more convenient to work in the Fourier space by defining Fourier transforms 
as, for example, 

+00 

C{k,oj)= j dt j dxe-^^''-'''^C{x,t). (4.12) 

— oo 

After grouping together all the terms that are linear in perturbation in the ex- 
pansion of equation (3.7), we obtain the linear stability equation 



(4.13) 



where 



C{k,u;,C) 



+ 



+ l-Ml(A- dolk^ + 



ga+A _ I ga-A _ i 

+ 1 + Ml] J2 e-(^+^/')(")^C(A;,t^) 



C{k,u;) 



m<n 



^ m>n 



-{M-l/l){m-n)\^m 



C{k,co). 



(4.14) 

Linear dispersion relations under general perturbations are solutions of (4.13), as 
we shall show below. 
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Before moving on, we reiterate that, although we have so far considered only 
straight steps, equation (3.7) also has steady state solutions corresponding to 
curved steps. In general, these solutions can be approximately evaluated, in 
the small amplitude limit, by expanding (3.7) to higher orders in (. The linear 
stability of these solutions can in turn be similarly investigated. However, these 

rgj 

calculations are in practice quite cumbersome. 
4.1. Stability of an Isolated Step 

As a special case, let us consider the linear stability of an isolated step. This 
situation is equivalent to the limit of infinite step spacing A = oo. Stability 
equation (4.13) simplifies to 

^ + 1 - (A - dolk^) lM{k, u) = (4.15) 
which gives the linear dispersion relation in the closed form 



= (A - d^lk^) \^fk'^ + /2 - 1 + (A - dQlk'^Y - 1 + (A - dQlk^Y (4.16) 



where the quantity I = \/l — A^/A is solved from (4.3). This dispersion relation 
is plotted in Fig. 2. 

It is evident from Fig.2 that, within certain parameter range, a straight 
step becomes linearly unstable against infinitesimal, long-wavelength fluctua- 
tions. This diffusive instability is of the same nature as that discovered by MuUins 
and Sekerka'^' for planar solidification fronts. The relevance of MuUins-Sekerka 
instability in step flow was flrst pointed out by Bales and Zangwill!^' The step 
problem in question here, although resembling closely that of solidiflcation, has 
important differences, in two respects. First, in solidiflcation, steady state planar 
fronts growing into supercooled liquid do not exist at general supercoolings, in 
contrast with the situation here where uniformly moving straight steps are valid 
solutions at arbitrary supersaturation. Second, a planar front in solidiflcation is 
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always linearly unstable to perturbations of sufficiently long wavelength, with the 
effect of surface tension merely stablizing the front on short lengthscales. But, 
again by contrast, straight steps can be completely stablized by strong enough 
line tension on all scales (see below) . Both differences between the two problems 
are due to the presence of the extra lengthscale, the adatom diffusion length Xg, 
in the step problem. 

It follows from (4.16) that, as long as do > a straight step is linearly 

stable against all infinitesimal perturbations, i.e., Re(ia;) < for all k. Since 
I = Vl — A^/A, this condition is equivalent to 

do > . (4.17) 

On the other hand, for smaller values of do, the straight step is unstable against 
long wavelength perturbations. The critical wavenumber below which perturba- 
tion grows is given by 



I 9 V2 

do . I do , d^ 



+ 



2A3/ V A3/ 4A6/2 



(4.18) 



To see if (4.17) can be satisfied in realistic experimental situations, we esti- 
mate various quantities using values typical to crystal growth from solution. As 
an example, usual crystal growth from solution has a value of supersaturation 
(coo -co)/co ~ 10~^. Further, we use ^Q/ksT = 25A, Ceq/cr = 10~'^,Xs = 800A, 
for NaCl crystals at T = 600K. This gives an estimate that do ~ 3 x 10~^, 
A fa 10""^, and therefore do/A'^ > 1. Thus (4.17) is most likely to hold for 
conventional crystal growth, implying complete linear stability of straight steps. 
However, instability may occur under MBE growth situations where much higher 
supersaturation can be achieved, as suggested in Ref.3. 

Although (4.17) is derived within the symmetric model, the qualitative con- 
clusion, that there exists a critical value of line tension above which the step is 
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linearly stable, is unaltered when asymmetric attachment kinetics are taken into 
account. Only the numerical value of the threshold (4.17) is modified. It will 
be shown later that in a one-sided model with strong attachement kinetics, the 
threshold (4.17) changes to do > A/2. So step instability is relatively easier to 
realize, even for the parameters of conventional crystal growth quoted above. 

4.2. Stability of Step Train Under Lateral Fluctuations 

We now examine the linear stability of an infinite train of equidistant steps. 
It is useful to distinguish two types of morphological perturbations: the lateral 
fiuctuations to the shapes of steps, and the longitudinal displacement of steps 
from steady state positions. Both kinds of perturbations are accounted for si- 
multaneously in the stability equation (4.13). 

In this section we first investigate the linear stability of steps with respect to 
fluctuations of the first type, assuming to this end that the average spacing be- 
tween steps is maintained at fixed value A. Still, shape perturbations on different 
steps need not be synchronous, but the most unstable situation corresponds to 
all steps having in-phase fluctuations. That is, (^{k,u>) = ('^(k^cj) for arbitrary 
m, n. In this case, equation (4.13) yields 



which can be solved implicitly for the dispersion relation u = iv{k). The presence 
of a zero mode in the spectrum, due to the translational invariance of the steady 
state solution, is evident by explicit inspection: a; = A; = solves (4.19). Un- 
fortunately, equation (4.19) does not afford an explicit analytic solution of (jj{k). 
However, one can still demonstrate that steps become completely linearly stable 




(4.19) 
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for strong enough line tension. The critical value of the dimensionless capillary 
length is determined via the condition 



k=0 



which gives the formula for the threshold line tension: 



do > 



a-\- — a.- 
(q!+ + a-Y 



A + 



Q+ A 



a_ A 



(ga+A _ 1)2 (ga_A _ i)2 



(4.20) 



where q;± = q;±(A, A) have been defined previously and can be solved from 
the steady state solution (4.3). Two limits, of large and small step spacing 
respectively, are analytically tractable. Specifically, we obtain from (4.3) that 



a+ = 1 + — A + 



(4.21) 



for A <^ 1 and 



1 + A 
1- A 



1 + 



1- A2 



/l-A X _ /l+A 



+ 



(4.22) 



for A » 1. Then condition (4.20) reduces, in the two limits, to 



(4.23) 



for A <^ 1 and 



A2/r~A^ AA 



(l-A)e" 



-(1 + A)e Vi-^'' 



+ 



(4.24) 



for A » 1. Hence the critical value of do needed to linearly stablize straight steps 
of finite spacing against lateral perturbations is always smaller than that required 
for an isolated step. This observation is consistent with the intuitive expectation 
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that smaller A suppresses lateral fluctuations of the steps, due to greater overlap 
of adatom diffusion fields. For asymmetric attachement kinetics, this conclusion 
remains qualitatively true as well, and has been reached by previous authors 
using quasi-static analyses. 

Our analysis is strictly linear, therefore does not address the inquiry regarding 
what happens once straight steps becomes linearly unstable. A natural expec- 
tation is that the steps may adopt profiles corresponding to other steady state 
solutions. A recent nonlinear, quasi-stationary analysis by Misbah and Rappel'"'' 
for the asymmetric model established the existence of cellular step solutions with 
a continuum band of possible wavelengths. Whether or which of these solutions 
can be actually selected is intimately related to their stability, an issue which 
remains unresolved. 

4.3. Stability of Step Train Under Longitudinal Fluctuations 

We now investigate the linear stability of an infinite step train against pertur- 
bations in the distances between steps. In this situation, an instability could lead 
to the coalescence or bunching of individual monoatomic steps into multiple-step 
bands, a feature frequently observed during crystal morphological changes. For 
simplicity, we assume that the line tension is large enough so that all steps remain 
flat so that we can set A; = in the stability equation. Arbitrary displacement 
of steps from their steady state positions can be decomposed into the normal 
modes as ('^{k^uj) = e^'^^'^Uq, with q being the real wavenumber. Without loss of 
generality, the wavenumber q can be further restricted to its flrst Brillouin zone 
q G j\- The stability equation (4.13) then reads 

+00 

An{iu)e"i^'' = Q (4.25) 

n=— oo 

where the coeffients Aniioj) are given by 

AqUu) = ^ + l- MfA + \ J ^ (4.26) 
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^n>o(^a;) = + 1 - ^) e-(^-i)-V' (4.27) 

A<o(^u;) = + 1 + ^) e(^+i)"^/' (4.28) 

where M{iuj) = lM{k = 0,uj) = Vl^ + 1 + iujl^. The frequency uj is in general 
complex, so it is appropriate to separate the dispersion relation into a real and an 
imaginary part, iuj{q) = ^R{q) +iO/((/). Positive values of imply exponential 
amplification of small perturbations hence linear instability. Once again one can 
check by explicit substitution the presence of a translational zero mode ^R{q = 
0) = ni(q=0) = 0. 

First let us examine in detail the analytically tractable limit of long- 
wavelength perturbations \qX\ ^ 1, corresponding to dispersion spectrum near 
the center of the Brillouin zone. Taylor expansion of the stability equation gives 
in this limit the relations 

^liQ) = -^^{qX) + OiqXf (4.29) 



where primes denote derivatives with respect to iu evaluated at iu = 0. Note 
that in the long wavelength limit, ^i{q) is linear but 0.ji{q) is quadratic in the 
wavenumber. Formula (4.29) can be evaluated to give 

(q;+ + Q!_)(q;+ sinh~^ + ct-sinh" 



4AP - Xa+ smh ^ ^ + Aq!_ smh ^ — 



The expression for the real mode is much more complicated and is not displayed 
here. Despite the complexity, behavior in the limits of large and small step 
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spacings can still be extracted, using the expansions (4.21) and (4.22). The 
results are, for A 1, 



nj{q) 



(1- A2)(l + A) 



(1-A2)(1 + A) 



-A2 



2A 





1 + A -./i±|A \ , 



(4.32) 



(4.33) 



and for A <C 1, 



niiq) = [A + OiX^)]qX 



(4.34) 



A2(1 + 3A2) 



A2 + 0(A4) {qX)\ 



(4.35) 



These asymptotic expansions are obtained from the full linear stability equation, 
hence involve no approximations. We note in particular that our result for fli{q) 



a step velocity depending upon only nearest terrace widths. Most importantly, 
however, in contrary to a pure neutral mode ^R{q) — found in their approxi- 
mate model, we obtain a weakly stable mode Qr < 0. 

While the magnitude of 0,ji{q) reflects the typical decay or amplification rate 
of perturbations, fi/(g) is related to the wave velocity for the propagation of 
disturbances. When the wavelength of disturbances is large compared with the 
step spacing, i.e., gA <^ 1, the group velocity of disturbances coincides with the 
phase velocity since ^i{q) is linear in q. Both are given by f = —Qj/q in the co- 
moving frame. We immediately see that for small step spacing X <^ 1, v = —V, 
so the disturbances remain stationary in the lab frame. On the other hand, for 
large step separation A » 1, v = 0, and the disturbances propagate along with 
the steps. 



agrees with that of Bennema and Gilmer,'"' obtained in a simple model assuming 
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The exact linear stability spectrum (4.25) over the complete Brillouin zone 
can only be solved niinicrically. However, we can attempt an approximate solu- 
tion utilizing a quasi-static approximation. As is usually done, this approximation 
assumes that the diffusion field adjusts quickly in response to the displacement 
of the steps. It amounts to neglecting the frequency dependence, or the memory 
effect, in the Green's function: M(iuj) = V^^ + 1 + iojp ~ VP + 1- Now (4.25) 
can be easily solved over the whole Brillouin zone to yield 



2 



A + 



1 



e^iri) v/ITT^cothAy^- 1 



(4.36) 



(1 — cos qX) 



ni{q) 



2VT+~PsmhX/l 



coth 



1 



+ 



1 



sin^A. (4.37) 



The quasi-static dispersion relation is shown, in dashed lines, along with the 
numerically determined exact spectrum, in Fig.3. The most important feature is 
that ^R{q) < for all q, suggesting complete linear stability of the equidistant 
step train against fluctuations in step spacings. 

To suppress the frequency dependence in M{iu>) when adopting the quasi- 
static approximation, we need |a;| <^ 1. From the quasi-static spectrum, it can be 
seen that typically Qj = 0(A) and 0,r = O(A^). Hence the condition of small 
supersaturation A <^ 1 is necessary. Another, more subtle, source of inaccuracy 
may occur in the quasi-static procedure, if the magnitudes of Qr and fi/ have 
large disparity between them. This seems to explain the discrepancies between 
the quasi-static and the exact spectra near the center of the Brillouin zone for 
small A in Fig.3. 
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5. GROWTH AND EVAPORATION WITH 
ASYMMETRIC ATTACHMENT KINETICS 



In the preceding discussions, steps are assumed to be ideal line sources/sinks 
near which a continuous adatom density is maintained at the local equilibrium 
value. This picture must be modified, in several respects, to accommodate to 
more realistic growth situations. First of all, effects such as the finite rate of 
adatom exchange at steps, insufficient concentration and non-ideal distribution 
of sink sites on steps, all lead to deviations of adatom density at the steps from 
the equilibrium value. In other words, growth rate is no longer solely con- 
trolled by adatom transport, but also by interface kinetics. Secondly, atomic 
exchanges on different sides of steps involve different kinetic barriers and are 
therefore nonequivalent. This asymmetry in attachment kinetics was pointed out 
by Schwoebel and Shipseyj^' and discussed extensively in subsequent literatures. 
Lastly, interface kinetics usually exhibit pronounced anisotropy. 

The effect of interface kinetics is usually taken into account by assuming an 
elevation of Cgtep from the equilibrium value linearly proportional to the adatom 
flux. Accordingly boundary condition (2.3) is replaced by 



where differences in the coefficients /?± measure the asymmetry of attachment 
kinetics. We have for simplicity neglected effect of anisotropy. In this section, we 
consider a situation with the most asymmetry, assuming local equilibrium near 
lower terraces but complete inhibition of adatom exchange from upper terraces. 
This amounts to setting /3+ = 0, /3_ = oo in (5.1). We further distinguish 
between two particular cases of crystal growth and evaporation, respectively, 
shown schematically in Fig.4. Again working in the moving frame with velocity 
V in the 2;-direction, the boundary conditions (3.2) and (3.3) are modified to the 




(5.1) 
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following dimensionless form 



(v + dx = -iVu)l ■ da""; (Vw)_ • do-" = (5.2) 



u+{x,C,t)^/\-dQK{C} (5.3) 

for the case of growth and 

{v + ^ j dx = {Vuy_ ■ dcr""] (V«)+ • do-" = (5.4) 



«_(x,C,t)=A-doAc{n (5.5) 

for the case of evaporation. The bulk diffusion equation (3.1) remains intact in 
both cases. Note however that for evaporation, the quantity A (> 0) is identified 
with the dimensionless undersaturation and the dimensionless diffusion field in 
this case is defined according to tt = (c — Coo)/cr in contrast with (2.4). The 
above models are analogies of the one-sided models in solidification. 

5.1. The Case of Growth 

Integro-differential equations for the step displacement are derived using the 
same method as before by applying Green's theorem to the domain shown in 
Fig.4. Since adatom density will not be continuous across a step, twice as many 
equations are needed compared to the symmetric case. A static version of the 
Green's function method was recently used by Misbah and Rappel'^°' in a bifur- 
cation study of the steady state solution of cellular steps. 

Consider an arbitrary point p e Ai on the terrace between steps A and B in 
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Fig. 4. Green's theorem gives 



u{p) = J dxiG{p\pi) 



A 



v{i-ui{pi)) + 



dU 



A 



+ j dcTiuiipi) ■ ViG{p\pi) 



+ J daiu^{pi)-ViG{p\pi) + V J dxiu^{pi)G{p\pi) 

B B 

(5.6) 

where to save notation the integration over dti is imphed in the above expression. 
When we let the point p approach step A, the second term on the r.h.s. with 
ViG develops an integrable singularity. After separating the contribution of this 
singularity, we obtain 

«i(C^,t) = 2 J dcnu%i)-ViG{p\pi)+2V J dxiu^{pi)G{p\pi)+ 



B 



B 



j dxi 

A 

X G{p\pi) 



dti 



dxi 



{x - xi))- 



t-tl 



(5.7) 



Another equation is similarly obtained letting p approach step B, 



u 



B 



(C^, t) = 2j damiipi) ■ ViG{p\pi) + 2 j dxi 



+ / dxiG{p\pi) 



^-(C^-Cf 



B 



dC[_ 
dxi 



{x - Xl)) 



V{1 - ui) + 



u-{pi) 



A 



dC(_ 
dti 



G{p\pi) 



t-ti 



(5.8) 



Once again, linearization of (5.7) and (5.8) leads to the stability equation. 
Here we first briefly mention the linear stability of an isolated straight step under 
lateral shape perturbations. In our formalism, it turns out the dispersion relation 
can again be solved closely, from the relation. 



iujf + 2 = (A - dolk^){lM + 1) + Al{kH + Vl + /2M - IM^) 



(5.9) 



where M{k,uj) = \Jk'^ + 1 + 1//^ + iu). The actual expression for a;(A;) is lengthy 
and not shown here. We only point out that the threshold line tension needed to 
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linearly stabilize the straight step is now given by 



do > 



2- A 



(5.10) 



while for smaller do the step is unstable to long-wavelength perturbations. Thus 
suppressing adatom attachment from the upper terrace makes the step more 
susceptible to shape fluctuations. Further, it can be shown that the effect of 
finite step spacing A is stablizing, which merely reduces the above threshold to 
smaller values. The calculation is straightforward and reproduces the results of 
previous authors.'^' 

Below, our major attention will be paid to the longitudinal modes of an 
infinite step train. For straight steps, the equations of motion are much simpler. 
Since = A and ({x, t) is independent of x, the integration over dxi can be 
carried out explicitly. And the equations for step A and B are: 



A = dh 



V{2 - A) + 2 



dti 



A 



+ 2 j (itiii:?(ti)[VG(C^,t|cf,ti)- VuG(C^,t|cf,ti) 



(5.11) 



u^_{t) 



dU 



A 



21^(1- A) + 2 



dt^ 



+ 2aJ dtiVuG{C'',t\Ct,ti) + V J dtiu^{ti)G{C'',t\C^,ti) 



(5.12) 



with the reduced Green's function 

g-(t-ti) 



{Z-Z1 + V{t - ti)f 



(5.13) 



Uniformly moving steady states are zeroth order solutions of (5.11) and (5.12). 
Setting = 0, (-^ = A, we obtain 



U- 



e-'^+^{a\ + A - 1) = e°-^(a2 + A - 1), 



(5.14) 



from which both the velocity of the step and the adatom density at the upper side 
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of the steps can be calculated. And the magnitude of discontinuity in adatom 
density across the step equals A — tt_. 

The linear stability equation, derived from (5.11) and (5.12) by linearization, 
follows the same structure as (4.25), but with only two terms n = 0, 1 in the 
summation. It reads 



M -1 
M + 1 



,{M-l)X/l 



iujf + 2 - A - 2M + AmVi + 



+ e 



-{M+l)X/l 



X 



iujf + 2 - A + 2M + AM^ 



= U- + a-l — 1) 



(5.15) 



where as before M{iu) = yl + W+iujP . It is easier in this case to extract the 
dispersion relation near the center of the Brillouin zone. We find, to leading 
order, that 



2A(2 - A)^ 



(1 - A)2(4-2A- A2) 

-A(2- A)3 



__2^A ^ 

e gA 



(1- A)2(4-2A- A2) 
in the limit of large step separation A ^ 1 and 



ni{q)^[A + 0{X^)] qX 



3Ai 



+ 0(A2) 



(5.16) 
(5.17) 

(5.18) 
(5.19) 
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in the limit of narrow step trains. The spectrum solved under the quasi-static 
approximation also takes a particularly simple form: 



^i{q) = T sinqX 



(5.20) 



2 p2\/l 



u_e 



-(1 — cos qX). 



(5.21) 



Note that both Qji and fi/ are of the the same order 0(A). In Fig. 5, a comparison 
of the exact spectrum with the quasi-static one suggests the latter to be a very 
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good approximation. Since i^R{q) < for all q, the infinite step strain is linearly 
stable against perturbations in step spacings. 

5.2. The Case of Evaporation 

With little modification, procedures used in the last subsection can be applied 
to the case of crystal evaporation. We will not repeat the analyses, but only quote 
the results. For an isolated straight step, the linear stability spectrum for lateral 
fiuctuations is the closed solution of 



iujl"^ + 2 = (A - dQlk^){lM - 1) - /^l{kH + y/lTl^M - IM^). (5.22) 



It follows that iuj is pure real and iu!{k) < 0, for arbitrary do and k. So a straight 
step is completely linear stable against shape fluctuations. This conclusion re- 
mains valid for equidistant step trains with finite A as well. 

The equations of motion for parallel steps resemble those in the case of growth 
and read 




(5.23) 



2 / dtiui{ti)[VG{C'',t\Ct.h)-VuG{C,t\Cf,h)l 




Bl 



G(C^,t|Cf,ti) 



(5.24) 



2A / dtiVi,G(C^,t|Cf,ti)-F / dti«^(ti)G(C^,t|Cf,ti). 



Stationary state adatom density and step velocity can be solved from 



u+ 



(1 + A - a' 



;2 ) = e"+^(l + A-al) 



(5.25) 
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and the linear stability equation is 



M + 1 



g(M+i)A/z r ^2 + 2 + A + 2M - AM vTT]2l + e-(^-i)V« 



X 



M- 1 



(5.26) 



+ 2 + A - 2M - AM^ 



«+e-*'^^M(M + a+/ + l). 



The quasi-static dispersion relation, 



2 -2X/1 



sin qX 



(5.27) 



A 




(1 — cos gA), 



(5.28) 



A 



is displayed in Fig. 6 with the exact spectrum. We immediately observe that 
^r{<i) > for all values of q. Therefore, while the infinite step train is linearly 
stable against lateral shape fluctuations, it is unstable to perturbations in step 
spacings, of all wavelengths. 



At finite temperatures, a monoatomic step on a singular crystal surface con- 
tains kink excitations. As the temperature is raised, the density of kink and other 
short-range structural excitations increases and the step can become molecularly 
rough. In fact, in our previous considerations of step kinetics, molecular rough- 
ness of steps is necessary for them to behave as good sources/sinks of adatoms 
so that the continuum description applies. Moreover, roughness of steps can per- 
sist to macroscopic scales comparable to step extension, since a step is always 
thermodynamically rough at finite temperature meaning that the amplitude of 
fluctuations diverges with the system size. Beside contribution from equilibrium 
fluctuations, step roughness is also affected by nonequilibrium properties of the 
growth process itself, for example, by the amplification of noises by the MuUins- 
Sekerka instability as we shall show below. 



6. MORPHOLOGICAL FLUCTUATIONS OF STEPS 
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Previous studies of step morphology were mostly confined to treatments of 

[12] 

the equilibrium fluctuations on weakly interacting steps. In addition, the effect 
of thermal and other excitations on the terraces are routinely neglected. It is not 
until recently, that attention has begun to be paid to step fluctuations under non- 

[13] 

equilibrium situations. In one such study by Uwaha and Saito, step roughness 
is examined by adding to the linear stability equation a hypothetical stochastic 
noise term representing the cumulative effect of all thermal fluctuations. The 
origin of their choice of the randomness, however, remains elusive. Another 
study, by Salditt and Spohn'^*' on the time-dependent step roughness, took into 
account bulk random noise from the terraces but neglected noise contributions 
from excitations on the steps. Neither investigation considers multiple steps. 

In this section, non-equilibrium step roughness is investigated using a 
Langevin formalism. Our approach is an extension of the general theoretical 
framework for the study of morphological fluctuations of solidiflcation fronts at 
or near equilibrium, developed by Langer'^^' and Karmaj^^' based on the proce- 
dures introduced by Cherepanova. In this formalism, non-equilibrium, hydro- 
dynamical fluctuations of the phase boundaries are evaluated as responses to the 
Langevin forces representing local thermodynamical fluctuations. We demon- 
strate that this Langevin formalism provides a consistent description of the ki- 
netic roughness of steps. The effect of finite step spacing on step roughness, and 
the correlation between fiuctuations on different steps, can all be conveniently 
evaluated. Our approach is limited in several aspects. The first and foremost lim- 
itation is related to the unresolved fundamental question of statistical mechanics 
of whether the Langevin formalism itself is applicable in treating statistical fluc- 
tuations of nonequilibrium systems. Second, our calculation is carried out in the 
small amplitude limit and allows neither overhangs nor crossings of steps. Sim- 
ilarly, no interactions between steps beside the pure diffusive coupling is taken 
into account. 

In the Langevin approach, stochastic noises are introduced to the diffusion 
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equation and to the Gibbs-Thompson boundary condition as 



^^ = V\-u + qM (6.1) 



u{x, C, t) = A- doK{C} - Pvn{C} + V{C} (6.2) 

where the two independent Langevin forces q and rj mimic fluctuations in adatom 
density on the terraces and near the steps respectively. The noises are assumed 
to be Gaussian distributed with zero mean and variances 

(^(r, t)qir', t')) = 2r(l - V')Sir - r')Sit - t') (6.3) 

{r]{x, t)r]{x', t')) = 2(3r5{x - x')5(t - t'). (6.4) 

Note that the two-dimensionaUy distributed noise q comprises two contributions: 
one coming from the randomness in the impinging flux of atoms from the vapor, 
and the other from the adatom terrace diffusion processes itseff. The magnitude 
r, as yet unspecifled, wifl be flxed by the requirement that the method recovers 
correctly the equilibrium fluctuations. 

Restricting ourselves to an inflnite equidistant step train in the symmetric 
model, we write the stochastic equation of motion for the nth step as 

A - dQK-l3vn + r/{r } - (J{x, C, t) 

= J2J dhdxiG {x, c\ t\xu cr, ^i) {v + ^) ^^-^^ 

where the two-dimensional noise is projected as 

(7{x, C, t) = J dtidriG{x, C, t\ri,ti)q{ri, h). (6.6) 

Considering small amplitude fluctuations around straight steps we linearize the 
above equation in the Fourier representation. The step proflles in response to 
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the noises satisfy 

+00 



J2 An-rn{k,io)C{k,io) = a^{k,u;)-f]^{k,io) (6.7) 

m=—oo 

where coefficients ^n-m only depend on the difference n — m: 

Aq = 1(a + 1 1 TTT I -^t^ ( 4t + /3 I -c^o/c^ (6.8) 



^n—m 



with M(A;,a;) as previously defined. Also, we have 



+00 



a^{k,u) = I d^ie-(^^-"^)/^-l^^-"^l^ ^^^^^^' 'j^ . (6.11) 

—00 

The infinite set of equations (6.7) have the Toeplitz form and can be solved 
straightforwardly using Fourier series. Defining 

CO 00 

A{e)= J2 ^(^)= (6-12) 



S{e) = J] ((7" - 77^^)6*'^^ (6.13) 
n=— 00 

we have 

A{e)z{e) = s{e). (6.14) 

Inverting the Fourier series, we obtain 

—■K 

which leads, noticing that {a'^a'^) only depends on n — m and that {rf'rf^) = 
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^nrm to the result 



-TT 



(6.16) 

This expression is the fluctuation spectrum for the infinite step train in response 
to stochastic noises from both the terraces and the steps. The stationary fluctu- 
ation spectrum can be evaluated from (6.16) by integrating out the arguments 
A;', uj' and uj: 

ere.) = / {r{ku)c\k'u')) . (6.17) 

To proceed further with (6.16), we shall adopt the quasi-static approximation 
as before by neglecting the uj dependence in M{k,(jj). Consequently, the noise 
correlators can be evaluated: 

{a''{kuj)a\k'J)) = T{27rfSik + k')Siu + a;')e-l-l^^i±^^^!^, (6.18) 

{rf{ku)rf{k'u')) = 2pr{2n)^6{k + k')5{uo + J). (6.19) 

After some algebra, we flnd the general expression for the stationary fluctuation 
spectrum of steps 

;o \ „ r / dee-'--'a{e) 



where 



a{e) = 2m+ a (6-21) 

^ ' ^ cosh MA -cos 6* ^ ' 



1 / 1 
h{e) =7 A - d^lk^ + 



I V e"+^ - 1 e"-^ - 1 Ml 

If l + Ml I -Ml 

+ 



I g(M+l//)A-i6' _ 1 ' ^{M-l/l)\+ie _ I 



(6.22) 
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0(9) = 1 + 2PM+ ^(^^,/,)\_,, _ ^ + ^^M-i/ll+iO _ 1 (6-23) 

and the quasi-static M = + k'^ + Equation (6.20) is the main result of 
this section. 

First it is important to examine what (6.20) suggests of the equihbrium fluc- 
tuations. Setting A = 0, / = cx), we have a{6) — c{6), b{6) — —dok^ and hence 
the equilibrium fluctuation spectrum 

4"C°.).. = J^. (6.24) 



eg 



We observe that the step roughness in equilibrium is independent of the step 
spacing A. And apparently, in equilibrium, fluctuations on different steps also 
decouple since (CaTC^^-A:) = for m 7^ 0. 



eg 

Under non-equilibrium growth situations, however, fluctuations on different 



steps become correlated since in general {^Ck^C-kj 7^ for m 0. Furthermore, 
equation (6.20) exhibits the familiar divergence of fluctuations at the onset of 
MuUins-Sekerka instability when b{9 = 0) = 0. To see this more clearly, we 
consider the simple limit of a single step corresponding to A = 00. In this case 
the kinetic roughness spectrum of the step reduces to 

^^^-^) " \dok^-A/l + l/Ml2\ ^^-^^^ 

where the denominator is proportional to the quasi-static dispersion relation 
of an isolated step. Extensions of the present calculation to include energetic 
interactions between steps and to incorporate asymmetric attachment kinetics 
may be warranted for direct comparison with experimental observations. 
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7. COLLECTIVE MOTION OF STEPS 
— THE PHASE FIELD METHOD 

A rich variety of collective behavior of steps are commonly encountered dur- 
ing changes of crystal surface topography. Steps are seldom uniformly spaced; 
pairwise grouping of steps, bunching of steps into macro-steps with height of mul- 
tiple atomic units, collision and annihilation of step with antisteps, have all been 
observed in various surfaces of, for example, alkali halides,'^*' semiconductors,'^^' 
and metals!^°' To address these effects, fully time-dependent treatments of the 
collective movement of steps are necessary. 

Theoretical analysis of the time-dependent multiple step kinetics turns out to 
be formidable. The most well known continuum approach is the kinematic wave 

[211 [221 

theory of Frank using the method of characteristics. Frank's theory treats 
bundles of steps as basic entities and monitors the temporal changes of the average 
step density. It does not trace the position of individual steps and is therefore 
more coarse-grained than the continuum BCF model. However, the method is 
found particularly useful in describing large scale changes of surface morphology 
involving macro-steps. Time evolution of step spacings in finite trains of steps 

[23] 

has been discussed by MuUins and Hirth, based on the assumption that the step 
velocity is a function of only two adjacent step intervals. Both approaches apply 
to straight steps only and do not allow lateral variation of step profiles. Other 
time-dependent studies '^"^ of step motion either follow the spirit of MuUins and 
Hirth, or resort to explicit microscopic Monte Carlo simulations of the growth 
processes. 

Meanwhile, the attempt towards a direct numerical solution of the original 
Stefan problem (2.5) to (2.7) is also difficult, hindered by the need to track 
explicitly all moving boundaries. This difficulty was highlighted in studies of 

[71 

dendritic solidification, and can be partly relieved using the so-called phase 
field method. It is the purpose of this section to develop a phase field model for 
the step problem. We shall demonstrate that this method provides a powerful 
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tool permitting detailed study of time-dependent step kinetics. 

The crux of the phase field approach lies on the introduction of an order 
parameter 0(r, t) indicating the phase at a particular position. In our model, 
local stable minima of the order parameter correspond to terraces whereas rapid 
spatial variation of the order parameter locates the position of steps. Now we 
introduce a stochastic phase field model in the following dimensionless form 



T^^^ = + a sin(7r0) -u + A + q^{r,t) (7.2) 

where the last equation can be derived through a pure relaxational dynamics 
Td(f)/dt — —SF/S(f) + from the hypothetical free energy 



F{(j),u} = J (fr 



t2 Q 

— (V0)^ + - cos(7r0) + (tt - A)0 

2 TT 



(7.3) 



The Langevin noise q satisfies (6.3) while has zero mean and variance 

(g^(r,t)g^(r',tO> = ^5(r - r')5(t - t')- (7-4) 

The mangitude of variance is so chosen as to reproduce the correct interface noise 
r] in the original sharp interface model discussed in the previous sections. 

Equations (7.1) and (7.2) are extensions of the classical phase field model 
for solidification to the present problem with thermal fiuctuations. A sinusoidal 
potential term is introduced to facilitate the description of multiple steps, by 
identifying the degenerate minima ~ (2i + l)7r ,where i is an arbitrary integer, 
with terraces. The location of the moving steps are defined by the condition 
(f){r,t) = 2z7r. The phase field model captures, phenomenologically, the effect 
of line tension and finite attachment kinetics by having finite values for ^ and 
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T, since ^ can be regarded as the thickness of steps and r reflects the rate of 
response of the phase field. Parameter a is the strength of the potential which 
can be taken as a constant of order unity, when considering the limit of small 
supersaturations. 

Formal methods have been used to establish the phase field model as a proper 
regularization of the original sharp interface problem (2.5) to (2.7). It has been 
shown,'^'' again formally using matched asymptotic expansions, that the phase 
field model recovers the sharp interface model by taking appropriate limits of 
quantities r ^ 0,^ ^ 0, and a -C 1. Following the standard procedures, it can 
be shown for our present model that the value of tt at a step satisfies 

Ustep ~ A - dQK - + T] (7.5) 

with the correspondences do ~ C'sA) ~ '^V^/Cj that the projected Langevin 
noise 77 is governed by the correlator (6.4). 

Equations (7.1) and (7.2), although stiff for small parameters ^ and r, are 

[27] 

suitable for direct numerical analysis. While simulations on the physical two 
dimensional geometry is straightforward, we focus here on the deterministic one 
dimensional problem, corresponding to straight steps. First of all, numerical 
evidence suggests that (7.1) and (7.2) support steady state solutions with con- 
stant velocity, for infinite equidistant step trains with arbitrary spacing A, and 
at all values of supersaturation < A < 1. The most interesting applica- 
tion, however, concerns the time-dependent motion of finite, in general non- 
equidistant step trains. Such a situation is illustrated in Fig. 7, where a snap- 
shot of the adatom density and order parameter profiles is displayed for a non- 
equidistant train of six steps. The parameters used in our numerical calculation 
are A = 0.1, a = 1,t = 2.5 x 10~^, ^ = 0.05, and the average spacing between 
steps is A 2.5. Strong overlap of the terrace diffusion fields is evident, even 
for average step spacing of A 2.5. To dispel a reasonable misgiving, it is not 
hard to realize that an effective short range repulsion exists between steps spaced 
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on the scale of ^. This conveniently avoids the possibility of steps behind over- 
taking others in the front, which in reality would create overhangs and involve 
prohibitive energy penalty. 

A train of finite number of steps can not propagate with equal spacing in- 
definitely, because there is no corresponding steady state solution. In Fig.8, 
we consider the temporal evolution of a train of six steps starting from a con- 
figuration of equal step spacing A = 0.5 at A = 0.1. Step configurations at 
t = 0.5, 20, 60, 100, 140 are shown in sequence. Both the leading step and the 
trailing step move faster than the steps inside the train, because each of them 
borders on an infinite terrace with richer adatom supply than other terraces. 
Therefore one expects the step train to pile up at the rear and spread at the 
front. This intuitive conclusion is indeed vividly observed in Fig.8, where we see 
the bunching of the last three steps into a triple step, and the continuous break- 
ing away of the leading steps from the train. Meanwhile, the overall extension 
of the step train increases monotonically in time. At t = 140, the total length of 
the train has increased from an initial value of 3.0 to 16.14, by a factor of five. 
It is also interesting to monitor the subsequent evolution of the triple step at 
the rear in Fig.8. First it is clear that a step with multiple height moves much 
slower than the isolated steps. Consequently, the triple step lags farther behind 
the other steps until the leading step within the bunch finds a larger terrace in 
the front and accelerates to break away from the group. After a further dissoci- 
ation of the remaining double step, the triple bunch eventually disassembles into 
elementary steps again. The qualitative conclusions we can now draw, that finite 
step trains are susceptible to bunching at the rear and spreading in the front, 
and that the average step density in a bunch decreases with time, are consistent 

[23] 

with the analysis of Mullins and Hirth and with the general results of Frank's 
kinematic wave theory !^^' 

In short we see that, the complex dynamics of a finite step train result from 
the combined action of two efi^ects — of step grouping from the rear and spread- 
ing from the front — each initiated from the two ends of the step train. Both 
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effects can nevertheless be more instructively illustrated when examined sepa- 
rately, in Fig. 9 and Fig. 10. For instance, Fig. 9 shows the temporal development 
of pairwise bunching undulations on a semi-infinite step train. Clearly, the con- 
figuration of the step train at any time consists of a region with disturbances 
and an unperturbed region, separated by a propagating "shock" front. Within 
the region of disturbances, a pairwise bunching of steps occurs. In addition, 
the front appears to propagate uniformly with a constant velocity. A different 
behavior is observed for a moving semi-infinite step train which terminates on 
the right, as shown in Fig. 10. Here no bunching instability is present and the 
step train spreads smoothly. It is then inferrable from the foregoing discussions, 
confirmed by our numerical simulations, that whenever a large step train of lower 
step density invades that of a higher density, a shock front froms at the interface 
which propagates into the train of higher density and leaves behind a region of 
bunched step configurations. For a high density train running into a low density 
train, on the other hand, step spacings at the interface vary smoothly in space 
and interpolate between values characteristic of the two step trains. 

These one dimensional examples which wc have hitherto considered are used 
to illustrate the utility of the phase field model to the study of microscopic ki- 
netics of monoatomic steps. The major limitation of the present method lies in 
its difiiculty to be generalized to the case of asymmetric step kinetics, due to the 
discontinuity in adatom density across the steps. Nevertheless, many physically 
important factors, such as the impurity obstruction to step motion, temporal 
changes of growth environment, spatial inhomogeneity in supersaturation, crys- 
talline anisotropy, and the effect of island nucleation, can all be conveniently 
investigated using this phase field approach. Detailed studies of these aspects, in 
two dimensions, will be presented in the future. 
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FIGURE CAPTIONS 



1) Schematic representation of a train of steps moving in the ^-direction. 

2) Linear stabihty spectrum (4.16) of an isolated step against shape perturba- 
tions, in the symmetric model. Curves correspond, from top down, to values 
of line tension (do) 0.25, 1, and 2 times the critical value A^Vl — A2/2. 

3) Dispersion relation of an infinite train of parallel steps against perturba- 
tions in step spacings, in the symmetric model. Solid lines are obtained 
numerically from (4.25), and dashed lines correspond to the quasi-static 
spectrum (4.36) and (4.37). Curves in the upper figure, from top down, are 
for A = 0.5, 1, and 2. Order of curves is reversed in the lower part. 

4) Schematic representation of step flow during: (a) crystal growth and, (b) 
evaporation with asymmetric attachement kinetics. 

5) Dispersion relation of an infinite train of parallel steps against perturba- 
tions in step spacings during crystal growth, in the one-sided model. Solid 
lines represent the numerically solved exact spectrum, and dashed lines 
correspond to the quasi-static spectrum. Curves, in increasing amplitude, 
are for A = 2, 1, and 0.5. 

6) Dispersion relation of an infinite train of parallel steps against perturbations 
in step spacings during evaporation, in the one-sided model. Numerically 
solved exact spectrum is shown in solid lines, and the quasi-static spectrum 
is shown in dashed lines. Curves, in increasing amplitude, correspond to 
A = 2, 1, and 0.5. 

7) Snapshot of: (a) adatom density and, (b) order parameter profiles, for a 
train of six steps. A = 0.1. 

8) Time evolution of a finite step train, starting from an initial configuration 
of six steps evenly spaced at A = 0.5. In the phase field (0) plot, time 
sequences are, from left to right, t = 0.5, 20, 60, 100 and 140. For clarity. 
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only three instances of the diffusion field (u) ai t — 0.5, 60, and 140 are 
shown. A = 0.1 

9) Pairwise bunching of steps at the rear of a semi-infinite step train with 
uniform initial spacing A = 0.75 at A = 0.2. Curves from left to right 
are at t = 24, 48, 72, 96, and 120. Adatom density profiles are vertically 
displaced by multiples of 0.2 to show clearly the propagation of bunching 
instability into the step train. 

10) Uniform spreading at the front of a semi-infinite step train with constant 
initial spacing A = 0.75 at A = 0.2. Step configurations from left to 
right correspond to t = 24, 48, 72, 96, and 120. Adatom density profiles are 
vertically displaced by multiples of 0.4 for clarity. 
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